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Abstract 

Perturbative unitarization from non-linear effects is thought to deplete the 
gluon density for transverse momenta below the saturation scale. Such ef- 
fects also modify the distribution of gluons produced in heavy-ion collisions 
in transverse impact parameter space. I discuss some of the consequences for 
the initial conditions for hydrodynamic models of heavy-ion collisions and 
for hard "tomographic" probes. Also, I stress the importance of realistic 
modelling of the fluctuations of the valence sources for the small- x fields in 
the impact parameter plane. Such models can now be combined with solu- 
tions of running-coupling Balitsky-Kovchegov evolution to obtain controlled 
predictions for initial conditions at the LHC. 
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Conditions 



1. Introduction 

The concept of perturbative saturation of the gluon density of a hadron at 
small light-cone momentum x was introduced originally to preserve unitarity 
of the scattering amplitude at high energy [l[ which is violated by linear 
perturbative QCD. In particular, non-linear processes such as "gluon recom- 
bination" should prevent a power-law divergence of the (unintegrated) gluon 
density <&(x, fcjj at small intrinsic transverse momentum k±, and instead lead 
to its "saturation": k±) ~ log l/kj_. 

McLerran and Venugopalan suggested that for a large nucleus and at 
small x that the problem of gluon saturation could be addressed by classical 
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methods [2J. In their model, the high density jj 2 of valence charges per unit 
transverse area on the light cone acts as a source for a classical Yang-Mills 
field of soft gluons. This classical field, indeed, was shown to exhibit satu- 
ration for kj_ < Q 2 ~ a 2 p 2 jsj]. Most important for the present purposes, 
however, is that perturbative gluon saturation also predicts a modified dis- 
tribution of small- x gluons in transverse impact parameter space: in the 
unitarity limit the local density of gluons, integrated over kj_, is no longer 
simply proportional to the density of valence sources as expected in the di- 
lute, linear regime. 

The emerging consequences for the initial conditions for hydrodynamic 
models of heavy-ion collisions were realized later 0,1^. Consider a collision of 
two heavy ions at non-zero impact parameter. Neglecting fluctuations of the 
local density of participants, their overlap area in the transverse plane has a 
short axis, parallel to the impact parameter, and a long axis perpendicular 
to it. If the produced gluons equilibrate then the pressure gradients convert 
this asymmetry of the initial density profile into a momentum asymmetry 
called "elliptic flow" , v 2 ~ (cos 20) . In the absence of any scales (such as the 
freeze-out temperature Tf, the phase transition temperature T c , or a non- 
vanishing mean free path A) hydrodynamics predicts that t>2 is proportional 
to the eccentricity e of the overlap area p, e = (y 2 —x 2 )/(y 2 +x 2 ). The 
average is taken with respect to the distribution of produced gluons in the 
transverse x-y plane. 

A simple initial condition assumes that by analogy to the Glauber model 
for soft particle production dN/dyd 2 r± ~ PpZt( r ±) = (Ppart( r ±)+Ppart( r -L))/ 2 > 
where p 1 t is the density of participants of nucleus i per unit transverse area. 
A ~ (5 — 20)% contribution of hard particles needs to be added in order to 
fit the centrality dependence of dN/dy; their transverse density scales like 
~ T A B(r±). 

High-density QCD (the "Color-Glass Condensate") predicts a different 
distribution of gluons in the transverse plane, corresponding to a higher ec- 
centricity e for intermediate impact parameters. In particular, when either 
A or B is dense the number of produced particles is proportional only to the 
density of the dilute collision partner, whose partons add up linearly. Hence, 
in the reaction plane, dN/dyd 2 rj_ ~ mm(Q 2 s A ,Q 2 sB ) ~ min(pp art , Pp art ) drops 
more rapidly towards the edge than dN/dyd 2 rj_ ~ Pp™ t [5(. Thus, a higher 
eccentricity is a generic effect due to a dense target or projectile. Specific 
numerical estimates at the limited available RHIC energy do depend on the 
model for the unintegrated gluon distribution, however. 
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2. Modelling fluctuations of the large-a? valence sources 



Numerical estimates for the density distribution of produced partons re- 
quire detailed modelling of the fluctuations of the large- a: sources in impact 
parameter space. For peripheral collisions, and for smaller nuclei such as Cu, 
this is obvious. However, differences of moments of the density distribution, 
such as the eccentricity e, exhibit sensitivity to fluctuations even for central 
collisions of heavy nuclei • This is due to the fact that in the presence of 
fluctuations e ^ for central b — y collisions (with x and y directions defined 
via the principal axes of the particle production zone), while it would oth- 
erwise vanish. The same applies to another moment, the "triangularity" jsj, 
which gives rise to ~ (cos 30) contributions to the azimuthal distribution of 
particles (a non-vanishing triangularity arises in the CGC framework away 
from midrapidity, due to evolution of Qf' B with y, even without fluctuations 
of the sources [5j). 

To model the fluctuations of the valence sources for the small-x fields 
in the transverse plane one notes that several distinct transverse distance 
scales are involved. The radius Ra of a nucleus is much larger than the 
radius R n of a nucleon (the confinement scale) and one may therefore treat 
their fluctuations classically; in other words, we consider a collision of two 
"bags" of nucleons (quite densely packed, though) instead of using a quantum 
mechanical wave function for the nucleus. The transverse coordinates of the 
nucleons can be sampled randomly from a Woods-Saxon distribution. Multi- 
particle correlations are usually neglected (see, however, ref. (of) except for a 
short-distance hard core repulsion which enforces a minimal distance ~ 0.4 fm 
between any two nucleons. 

The scale where particle production occurs within the CGC framework is 
~ 1/Qs- It is again much smaller than the radius of a nucleon, 

^- < R N < R A , (1) 

and so the fluctuations of the valence charge density within a nucleon could 
again be treated independently from particle production. However, in prac- 
tice nucleons are usually treated as hard spheres with uniform density. 

Once a configuration of valence charges in the transverse plane has been 
obtained by Monte-Carlo methods, one constructs the unintegrated gluon 
distribution (the small- a; fields) $(x, k\\ r±) at every point r_|_ in the trans- 
verse plane [10|. The large- £ valence charges can still be treated as frozen 
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sources since their fluctuations in the transverse plane occur over time scales 
much larger than the small- a; evolution. In the model of Kharzeev, Levin 



and Nardi [111 , for example, the unintegrated "Weizsacker- Williams" gluon 



density of a nucleus (per unit transverse area) is given by 

k h r±) - ® { ?'' r± \ , 2 , • (2) 

a s max(Qi(x] r±), k z ± ) 

Ql(x; r±) = Ql(r±) (x /x) x exhibits the growth of the saturation momentum 
with energy expected from quantum evolution. In turn, Qg(r^) corresponds 
to the density of valence charge (squared) at the scale x and needs to be 
determined for each configuration individually. 

The fluctuations in the distribution of the hard sources in the trans- 
verse plane should not be confused with their distribution in color space. In 
the CGC framework, the color charge density per unit area of a nucleus or 
hadron is a stochastic variable, and all observables need to be averaged over 
its distribution. In the MV model, for example, this distribution is a local 
Gaussian, 

^]~exp(-jdW^^) . (3) 

Since there are infinitely many points in a transverse domain on the order of 
a fluid cell (linear dimension ~ ^/T) these MV "fluctuations" of the valence 
charge density do not correspond to fluctuations of the initial condition (par- 
ticle density) for hydrodynamics. While quantum evolution in x modifies the 
MV weight functional to a non-local distribution, the correlation length nev- 
ertheless remains of order 1/Q S and therefore much smaller than the size of a 
fluid cell. Fluctuations of the evolution ladders should therefore have a small 
effect on the hydro dynamical evolution and on "global" event properties (v 2 
etc.). 



3. Applications 

Fig. [T] shows the elliptic flow v% measured in heavy-ion collisions at RHIC 



scaled by the eccentricity e of the overlap zone [12|. As already mentioned 
above, in the absence of any scales such as a non-zero mean free path, v^/s 
would be independent of the transverse density of particles. Indeed, if the 
t> 2 data is scaled by the eccentricity obtained from a CGC implementation 
based on the KLN model then the required breaking of scale invariance is 
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Figure 1: Left: t^/e versus the transverse density [12|; i>2/£cgc nas been scaled by 1/2 
for better visibility. Right: Correlation between the suppression of hard non-photonic 
electrons and v 2 of bulk particles expected from a gravity dual model [17| ■ 



seen to be lower than for Glauber-like initial conditions. Actual solutions of 
viscous hydrodynamics (for V2) appear to confirm this simple observation in 
that the slope of v 2 /e versus transverse density is sensitive to the distribution 
of produced particles (lif . State-of-the-art simulations with realistic QCD 
equation of state and with initial conditions which account for the above- 
mentioned fluctuations will provide further insight. 

Other recent studies jsl, H3] have shown that fluctuations in the initial 
state which evolve through the hydrodynamic expansion can give rise to 
various structures observed in two-particle angular correlations, such as the 
"near-side ridge" or the "away-side cone". The fact that these correlations 
are very l ong range in rapidity provides a link to CGC physics of heavy-ion 



collisions [15|: the classical color fields introduced with the MV model are 
boost invariant, while quantum evolution in logl/x diminishes correlations 
only over rapidity intervals exceeding ~ l/a s . 

The modification of the distribution of produced particles in r±_ space in 
the non-linear regime affects not only the hydrodynamical evolution but also 
energy loss of hard probes. This observable provides a "tomographic probe" 
of the density distribution of the bulk. The higher eccentricity of CGC-like 



initial conditions increases the azimuthal asymmetry of high-p^ jets 16]. On 



the other hand, at fixed multiplicity dN / dy and v% of bulk particles it reduces 
the suppression of hard non-photonic electrons expected from a gravity dual 
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model, c.f. fig. Q] (right) [13]. 



4. Future developments 

Present CGC based estimates of the initial conditions for hydrodynamics 
from fluctuating valence sources usually rely on simple KLN-like models for 
the x, k± and r±_ dependence of the unintegrated gluon distributions (uGD). 
To improve the reliability of the approach one should incorporate constraints 
resulting from the measured (centrality dependent) ^-distributions in d+Au 
collisions at central and at forward rapidity. Furthermore, to actually test the 
theory of small- £ evolution it is important that extrapolations to LHC energy 
are based on uGDs which truly solve those evolution equations. Significant 
progress has been made recently in solving the running-coupling Balitsky- 



Kovchegov equation [18j, and those results can now be combined with the 



MC codes which generate configurations of large- x valence sources [10 . 
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